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Abstract 

We study the production of and Z gauge bosons in proton-proton 
and proton-neutron collisions up to a center-of-mass energy y/s = 0.5 TeV, 
in connection with the realistic possibility of having proton-deuteron colli- 
sions at RHIC with a high luminosity. We stress the importance of measuring 
unpolarized cross sections for a better determination of the flavor asymmetry 
of the sea quarks. Since polarized proton beams will be available at RHIC, 
we evaluate several spin-dependent observables like, double helicity asym- 
metries All as a test of the sea quark polarization, double transverse spin 
asymmetries Att as a practical way to determine a new structure function, 
the nucleon transversity distribution h±(x), and parity- violating asymmetries 

as a further test of our present knowledge of polarized parton distribu- 
tions. 
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1. INTRODUCTION 



A Relativistic Heavy Ion Collider (RHIC) is now under construction at 
Brookhaven National Laboratory and, already more than two years ago, it 
was realized that one should propose a very challenging physics programme M , 
provided this machine could be ever used as a polarized pp collider. Of cour- 
se all these considerations relie on the foreseen key parameters of this new 
facility, i.e. a luminosity of 2.10 32 cm~ 2 sec~ 1 and an energy of 50 — 250 GeV 
per beam with a polarization of about 70%. Since then, the RHIC Spin 
Collaboration (RSC) has produced a letter of intent ^ and has undertaken 
several serious studies in various areas which have led to a proposal^ which 
has now been fully approved. 

Among the very many basic hadronic reactions which will be studied in 
this experimental programme, we recall that jet production and direct pho- 
ton production are sensitive to the size of the gluon polarization, whereas 
Drell-Yan lepton pair production will allow to determine the magnitude and 
the sign of the sea quark polarization^ 1 ' 4 !. Recent measurements of proton 
and neutron Fi structure functions by the NMC at CERNt 5 ! have yielded 
a violation of the Gottfried sum rule^ 6 ' whose most natural interpretation is 
the evidence for a flavor asymmetry in the light sea quarks (i.e. d(x) > u(x)). 
In a recent paper! 7 !, we have shown that this can be studied further in the 
production of W ± and Z bosons in a proton-proton collider because it is 
dominated by quark-antiquark annihilation and therefore sensitive to the sea 
distributions. It was stressed that in the central region the predicted unpo- 
larized cross sections can vary by a large factor depending on, whether or 
not, one assumes a flavor symmetric sea. The measurement of charge asym- 
metries in the W + and W~ production was also suggested as an interesting 
possibility^ 8 '. Parity violating asymmetries with either one beam polarized 
Al or with two beams polarized have been calculated I 7 ! and we showed 
that their measurement, will give us a good calibration of the polarization 
of quarks and antiquarks, for both u and d. In this paper we are further 
investigating and Z bosons production and, in particular, we emphasize 
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some new possibilities if one can measure simultaneously proton-proton and 
proton-neutron collisions, assuming the feasibility of tagging the spectator in 
proton-deuteron interactions ™. We will also consider several spin dependent 
observables, with longitudinally and transversely polarized proton beams, 
which contain useful information on the polarized parton distributions. 

The paper is organized as follows. In section 2 we compare unpolarized 
cross sections in pp and pn collisions and we show it allows a much improved 
determination of flavor asymmetry of the sea quarks. In section 3 we evaluate 
double helicity asymmetries All for pp collisions with both proton beams 
polarized and parity violating asymmetries Al for pn collisions with only one 
proton beam longitudinally polarized. Finally in section 4 we consider the 
situation of pp collisions with both proton beams transversely polarized and 
we will see that only in the case of Z production, one can extract the new 
structure function h\(x) ^ 10 \ the transversity distribution, which indeed mea- 
sures the correlation between the left-handed and the right-handed quarks 
in a transversely polarized proton. We will give our concluding remarks in 
section 5. 

2. UNPOLARIZED CROSS SECTIONS 

Let us recall that the differential cross section for the reaction 

pp — > W ± + anything (1) 

can be computed directly in the Drell-Yan picture in terms of the dominant 
quark-antiquark fusion reactions ud — > W + and ud W~ and we have for 
the Standard Model W + production 

= G F nV2r ■ - [u(x a , M^)d(x b , M&) + (u <-> d)] (2) 

with 

G f = 7^77^—. i n ' x & = Vre y , x b = ^/re~ y , and r = M^/s . 
V2M^ sin 2 6 W 

Here sin 2 9w is the weak mixing angle and quark flavors are interchanged 
in eq.(2) for W~ production. Clearly these y distributions are symmetric 
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under y — > —y. For the reaction 

pn — > VF =I= + anything (3) 

under the assumption of isospin invariance for the nucleon, i.e. u p (x) = 
d n (x) = u(x), etc... (also for antiquarks), one has for W + production 

-^f- = G F nV2r ■ - [u(x a , M^)u(x b , M 2 ,) + d(x a , M^)d(x b , M 2 ,)] (4) 
ay 3 

which is no longer symmetric under y — > — y and it simply follows that for 

W~ production one has 

-Hy- {y) = ^- { - y) - (5) 
All these cross sections will strongly be dependent on the antiquark 
distributions whose present determination is more uncertain than that of 
u{x) and d{x). Following ref.[7] we will take the valence quarks u v {x,Q 2 ) 
and d v {x, Q 2 ) from ref.[ll] with the appropriate evolution from Q 2 = 4 GeV 2 
to Q 2 = M^y. We will assume that the antiquark and sea quark distributions 
go like (1 — x) 10 , flavor symmetry breaking being obtained by taking 

x5q(x) = xu(x) — xd(x) = ax b (l — x) 10 (6) 

with a = —0.22 and b = 0.45. For our calculations we have used M\y = 

80.22 GeV, the latest value given in the Particle Properties Data Booklet. 

This value is lower than the one used in ref. [7] and it leads to slightly larger 

predictions for pp collisions. The results are shown in figs. 1 and 2 for two 

different center of mass energies y/s = 350 and 500 GeV, the lower one being 

the maximum energy in pn collisions when the proton and the deuteron beam 

momenta are 250 GeV jc. Given the luminosity mentioned above, we see that 

one expects a fairly large event rate for all these cases. As in ref. [7] we could 

have studied the sensitivity of these cross sections to various choices of the 

antiquark distributions, but in order to avoid some uncertainties related to 

their absolute normalization and to pin down the flavor asymmetry of the 

sea we introduce the following quantity 

= daf p + jdy - daf n + /dy + daf p ~/dy - daj^/dy 
W ~ d*% + /dy + da^ n + /dy + da^~ jdy + da^jdy ' 1 J 
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In terms of parton distributions, it simply reads 

5q(x a )5q(x b ) + (x a <-> Xb) , . 

rtw = 7 x — x 7 \ [o) 

aq{x a )aq(x b ) + (x a <-> x b ) 

where ^ j (/(x) = it (a;) ± and similarly for antiquarks. Clearly R\y is 
symmetric under y — > — y and i?vK = if the sea is flavor symmetric, i.e. 
u(x) = d(x) that is Sq(x) = 0. The result of the calculation for our choice of 
the flavor symmetry breaking (see eq.(6)) is shown in fig. 3 for two different 
energies and it provides a very clear test for the antiquark distributions. We 
find that Rw is positive which is due to the fact that for our choice Sq(x) < 0, 
because d(x) > u(x). 

In the case of Z production, similar calculations can be done in the 
Drell-Yan picture in terms of uu and dd annihilations. For pp collisions we 
have 



= -^f- -3^^+ b ^ M z)Qifrb, M 2 Z ) + (x a <- x b )] (9) 



d °PP 
dy 

i=u,a 



where and bi denote the vector and axial vector coupling constants of the 
quark i to the Z boson. For pn collisions one obtains the corresponding 
expression da^jdy by making the appropriate substitutions as above for the 
case of production. The results are shown in figs. 4 and 5 for Mz = 
91.173 GeV at two different energies where we see, as expected, a symmetric 
distribution for pp and asymmetric for pn collisions. 

Let us also consider a quantity which is less sensitive to absolute nor- 
malization uncertainties i.e. 

Rz = dcr^ p /dy - da^Jdy 
Z da^ p /dy + daj^/dy 

and which reads in terms of the parton distributions taken at Q 2 = M§ 

Rz ={ (al + bl) [u(x a )Sq(x b ) + u(x a )5q(x b )] 

-(a 2 d + b 2 d ) [d(x a )5q(x b ) + d(x a )5q(x b )] } / 

( a i +b ^ [qi(x a )vq(xb) +q i (x a )aq(x b )} (11) 

i=u,d 



In this case, even for Sq = we don't have Rz = 0, because the couplings 
of u and d quarks are different and we find Rz negative because [a 2 , + b 2 ^ > 
(tty + b^) ■ The results are given in figs. 6a and b for two different energies 
where the dashed curves correspond to 5q = and solid curves to our choice 
of the flavor symmetry breaking (see eq. (6)). Although the Z production 
will be less copious this comparison between pp and pn collisions will also 
allow to discriminate between Sq = and 8q ^ 0. 

3. HELICITY ASYMMETRIES WITH LONGITUDINALLY 
POLARIZED PROTON BEAMS 

Since RHIC will be built to be used as a polarized pp collider, let us 
now investigate what we can learn from the measurement of various spin- 
dependent observables. We first start with one longitudinally polarized pro- 
ton beam and we consider the parity-violating helicity asymmetry defined 
as 

Here cr^ denotes the cross section where the initial proton has helicity h. 
This asymmetry will be expressed in terms of the parton helicity asymme- 
tries Af(x, Q 2 ) = f+(x, Q 2 ) — f-(x, Q 2 ) where f±(x, Q 2 ) denote the parton 
distributions in a polarized nucleon either with helicity parallel (+) or an- 
tiparallel (— ) to the parent nucleon helicity. In the Standard Model the W 
is a purely left handed current and this asymmetry reads simply, for W + 

production in pp collisions, 

Au (xg, M^)d (xb, Mgj - (u 
L[V) u(x a ,M^)d(x b ,M^) + (u^d) 
assuming the proton a is polarized. For W~ production quark flavors are 

interchanged. This asymmetry, together with two other parity-violating he- 
licity asymmetries related to it, have been studied in ref. [7], so here we will 
rather consider the case of pn collisions which will be also combined with 
some of our previous results on pp collisions from ref. [7]. For W + produc- 
tion in pn collisions one has 

= Au (x a , Mgj u (x b , Mgj - Ad (x a , M$y) d (x b , Mgj 
u (x a , 
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assuming the proton is polarized* 1 ). In order to calculate these asymmetries 
we have to use a model for the various partons (valence quark, sea quark, 
antiquark) helicity asymmetries Au v , Au 8 , Aw, Ad v , Ad 8 , Ad compatible 
with polarized deep inelastic scattering data. We recall that in leptopro- 
duction with both lepton and nucleon longitudinally polarized, one extracts 
the structure function gi(x) = 1/2 [Aq^x) + Aq^x)] which is, so far, 

i 

our unique source of information on the quark and antiquark helicity asym- 
metries. Many parametrizations have been proposed in the litterature, but 
following ref. [7] we will adopt a recent suggestion! 12 ] to construct polarized 
quark distributions from unpolarized quark distributions, based on the exten- 
sive use of the Pauli exclusion principle. In this approach^ 2 ) one successfully 
relates the violation of the Gottfried sum rule to the Ellis- Jaffe sum 

ru l e [i4] 

defect reported earlier by the EMC in the proton case! 15 ! . So concerning the 
valence distributions we will take 

Au v (x) = u v (x) — d v (x) and Ad v (x) = —1/4 d v (x), (15) 

whereas for sea quark, and antiquark, assuming that Aq s (x) = Aq s (x), for 
the H's we will take following ref. [7] 

Au(x) = u(x) — d(x) = 5q(x) (16a) 

which leads to a large negative u polarization. Concerning the d's, in ref. [7] 
their polarization was taken to be small and positive, but our analysis^ of 
the recent SLAC data on polarized neutron deep inelastic scattering! 16 ] i s 
consistent with 

Ad(x) = 0, (166) 

and this will be assumed from now on. We show in fig. 7 our choice of the 
various parton polarizations Aq/q taken at Q 2 = M^. Returning now to Al, 
we show the results of our calculations for pn collisions in figs 8a, b at two 

In the proton-deuteron collisions we assume that only the proton beam is polarized, 
because it is not clear to us that it will be possible to have polarized deuteron beams. 
( 2 ) For a more recent work on the description of quark distributions in terms of Fermi- 
Dirac distributions see also ref. [13]. 
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different energies and we observe that there is very little energy dependence. 
The general trend of A L is similar to what we obtained in ref. [7] for pp 
collisions. It can be understood from eq. (14) and we see that for y = —1, Al 
is sensitive to the antiquark polarizations since 

A?*--** and AT-~^ (17) 
d u 

whereas for y = +1, it is sensitive to the quark polarizations since 

and AY~~¥- (18) 
L u L d v ' 

So it is easy from fig. 7 to anticipate the results shown in figs. 8a, b and 
in particular if Ad is not equal to zero, it will clearly show up in A L V+ 
near y = — 1. From the measurement of the cross section da/dy and the 
parity-violating asymmetry Al for the four cases pp — > and pn — > 
one can obtain a quantity which is less sensitive to absolute normalization 
uncertainties, that is 



a L 



„ (A L da/dy)j; - (A L da/dy)X + (A L da/dy)^ p - - {A L da/dy)^ n 



daW + /dy + daW + /dy + d*%~ /dy + daf n ~ /dy 

(19) 

In terms of parton distributions taken at Q 2 = it simply reads 
w _ aAq(x a )5q(x b ) + aAq(x a )5q(x b ) 

(It — ": r — — ; r z r v^LH 

aq{x a )aq{x b ) + [x a <-> x b ) 

where aAq(x) = Au(x) + Ad(x) and similarly for antiquarks. So in a situ- 
ation where the sea is flavor symmetric i.e. 5q(x) = and the antiquarks 
are unpolarized i.e. Aq(x) = 0, one has = 0. Our results are shown in 
figs. 9a, b for two different energies. The solid curve corresponds to 8q(x) ^ 
and Aq(x) ^ so is positive and lies between 5% and 15%. However if we 
assume flavor symmetry breaking for the sea i.e. Sq(x) ^ 0, but unpolarized 
antiquarks Aq(x) = 0, we get the small dashed curve, which compared to 
the solid curve exhibits clearly the sensitivity to Aq(x) near y = —1. Finally 
one can also consider the case where 5q(x) = so Au(x) = and Ad(x) ^ 
at variance with eq. (16b), then we would have ^ and of opposite sign 
to Ad(x). 



We now turn to the case of Z production where similar calculations can 
be done. For A L we obtain the following expression 

A _ [2a u b u Au(x a )d(x b ) + 2a d b d Ad(x a )u(x b )] - (q <-> q) 
[(a 2 u + b 2 u )u(x a )d(x b ) + (a 2 + b 2 )d(x a )u(x b )]+(q^q) 

The results are shown in fig. 10 at two different energies. We get large positive 
values near y = +1, whereas it is almost zero near y = — 1, at variance with 
the case of pp collisions where it was also large and positive near y = — 1^ 
and this is related to the fact that d(x) > u{x) as we will see now. Like for 
W ± production combining pp and pn measurements, we now consider the 
following quantity 

a z = ( A Lda/dy)p p - (A L da / 'dy) z pn 
ttL " ' da^/dy + daZJdy 

In terms of the parton distributions taken at Q 2 = Mf it simply reads 

z _ [2a u b u Au(x a ) - 2a d b d Ad(x a )] [u(x b ) - d(x b )] - (q^q) 
^ ~ M + K) u{x a ) + (a 2 + bj) d(x a )] [u(x b ) + d(x b )} +(q~q)' 

Clearly in a situation where the sea is flavor symmetric and the antiquarks 
are unpolarized one has a£ = 0. We show our results in figs 11a, b for two 
different energies. The solid curve corresponds to Sq(x) ^ and Aq(x) ^ so 
af is around +10% near y = — 1 and reflects the trend of A L in pp collisions 
whereas af is around -10% near y = +1 and is driven by Al in pn collisions. 
For unpolarized antiquarks Aq(x) = 0, we get the small dashed curve which 
goes to zero near y = — 1 as Al in pp collisions for this case, but it follows 
the solid curve near y = +1. Finally we can also consider the case where 
Sq(x) = 0, so Au(x) = and Ad(x) ^ ; then we would have af ^ and 
with the sign of Ad(x). 

Before closing this section let us consider, in pp collisions where both 
proton beams are polarized, another observable which is very sensitive to 
antiquark polarizations, that is the parity-conserving double helicity asym- 
metry All defined as 

A = + do-— - aV + - - aV-+ _ 
da ++ + da-- + da^ + da |_ 
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Here crh 1 h 2 * s the cross section where the initial protons have helicities hi 
and hi- This asymmetry, in W + production, reads simply 

Au (x a , Mjj Ad (xb, M^) + (m ~ rf) 

For production quark flavors are interchanged. It is clear that All(v) = 
All{—v) and that All = if the antiquarks are not polarized, i.e. Au(x) = 
Ad(x) = 0. Similarly for Z production we find 



K + 6?) [Aft (x a , Mj) Aq, (x b , M|) + (x a <-> x 6 )] 



i=u,d 



J2 (a? + 6?) [ft {x a , M 2 Z ) q z (x b , Mf) + (x a <-> x b )] 

i=u,d 

(26) 



which will vanish for unpolarized antiquarks. We show in figs 12a, b our 
predictions for the three cases at two different energies and we now try to 
understand these results. Clearly as a consequence of eq. (16b) A L l = for 
W + production but if Ad(x) ^ 0, it would be non-zero and of opposite sign 
to Ad(x). For W~ production for y = we get All = —^f^r evaluated at 
x = M w /y/s which gives around +10% according to fig. 7. From the trend 
of the d and u polarizations shown in fig. 7 we also expect All to be almost 
constant for — 1 < y < +1. Finally for Z production as a consequence of 
eq. (16b) All does not depend on the d quark polarization. For y ~ ±1 
one has approximately All ~ - A ^ x ^ evaluated at x a ~ 0.65 and 

x b ~ 0.10 for y/s = 350 GeV that is, according to fig. 7, All h40%. 

For y/s = 500 GeV we get All ~ +25%. The values at y = are smaller 
because of the maximum of the unpolarized cross section. 

4. DOUBLE SPIN TRANSVERSE ASYMMETRY WITH 
POLARIZED PROTON BEAMS 

So far we have considered collisions involving only longitudinally polari- 
zed proton beams, but of course at RHIC transversely polarized protons will 
be available as well^. This new possibility is extremely appealing because 
of recent progress in understanding transverse spin effects in QCD, both at 



leading twist! 10 ! and higher twist levels! 17 ^. For the case of the nucleon's 
helicity, its distribution among the various quarks and antiquarks can be 
obtained in polarized deep inelastic scattering from the measurement of the 
structure function gi(x) mentioned above. However this is not possible for 
the transversity distribution h\{x) which describes the state of a quark (an- 
tiquark) in a transversely polarized nucleon. The reason is that hi(x), which 
measures the correlation between right-handed and left-handed quarks, de- 
couples from deep inelastic scattering. Indeed like gi(x), h\(x) is leading - 
twist and it can be measured in Drell-Yan lepton-pair production with both 
initial proton beams transversely polarized^ 10 !. Other possibilities have been 
suggested^ 18 ! but in the framework of this paper, we will propose a practical 
way to determine h\{x) using gauge boson production in pp collisions with 
protons transversely polarized. Let us consider the double spin transverse 
asymmetry defined as 

da n - da n 

a tt = —j {ft) 

da + da 11 

where o"|f (^TJ.) denotes the cross section with the two initial protons trans- 
versely polarized in the same (opposite) direction. Assuming that the under- 
lying parton subprocess is quark-antiquark annihilation, we easily find for Z 
production 



{bi ~ a 2 t ) Ihfixa^ixb) + (x a <-> x b ) 



, i=u,d 

Wt = 



(28) 



{ a 1 + h 1) [Qi(xa)qiixb) + (x a <-> x h )\ 

i=u,d 

This result generalizes the case of lepton-pair production! 10 ] through an off 
shell photon 7* and corresponding to bi = and = e^, the electric charge 
of qi. For W ± production, which is pure left-handed and therefore does not 
allow right- left interference, we expect Att = 0, since in this case of = b\. 
This result is worth checking experimentally. 

We recall that h\{x) has never been measured, so to evaluate Att one 
has to make an assumption on the possible magnitude of h\(x) for quarks 
and antiquarks. According to the MIT bag model! 10 ! one expects relativistic 
effects to increase the magnitude of hi(x) and since in the non relativistic 
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quark model one has 

h\(x) = Aq(x) (29) 

we will speculate that eq. (29) holds for quarks and antiquarks. The results 
of our calculations are shown in fig. 13 at two different energies and we see 
that Att has a similar trend as All in figs. 12a, b, as a consequence of 
eq. (29). Clearly this prediction is only a guide for a future experiment at 
RHIC which will indeed lead to the actual determination of hi(x). 

5. CONCLUDING REMARKS 

We have seen that W ± and Z production in pp and pn collisions at 
RHIC up to y/s = 500 GeV will provide an extremely valuable source of 
information on parton distributions. The high luminosity of the machine will 
allow copious production of these gauge bosons and the precise measurement 
of unpolarized cross sections will answer unambigously the relevant question 
of flavor symmetry breaking of the light sea quarks (or antiquarks) of the 
nucleon. Since polarized proton beams will be available at RHIC with a 
high polarization, it will be possible to use longitudinally polarized beams 
to measure parity-violating asymmetries Al- These asymmetries are large 
and contain clear signatures of the properties (magnitude and sign) of parton 
polarizations and it will greatly improve the calibration of light quark and 
antiquark polarized distributions. The absence of antiquark polarizations 
implies a vanishing double helicity asymmetry All, which is rather easy 
to test. Of course, all these informations will be complementary to the data 
extracted from polarized deep inelastic scattering experiments both on proton 
and neutron targets and will supply a detailed knowledge of the quark spin 
structure of the nucleon. Finally the use of transversely polarized proton 
beams for Z production in pp collisions will allow the first determination of 
the transversity distribution h\{x) for quarks and antiquarks, which is so 
far totally unknown. To conclude, we hope the arguments given above are 
strong enough to make gauge boson production a vital part of the future 
experimental programme at RHIC with polarized proton beams. 
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FIGURE CAPTIONS 

Fig. la da/dy versus y for W + and W~ production in pp collisions at 

= 350 GeV. 
Fig. lb Same as (a) at y/s = 500 GeV. 

Fig. 2 da/dy versus y for W + production in pn collisions at y/s = 350 
and 500 GeV. W~ cross section is obtained by symmetry around 
y = 0. 

Fig. 3 The ratio Rw (see eq. (7)) versus y at y/s = 350 and 500 GeV. 
Fig. 4 da/dy versus y for Z production in pp collisions at y/s = 350 and 
500 GeV. 

Fig. 5 da/dy versus y for Z production in pn collisions at y/s = 350 and 
500 GeV. 

Fig. 6a The ratio Rz (see eq. (10)) versus y at y/s = 350 GeV. Dashed 
curve corresponds to 5q = 0, solid curve corresponds to the choice 
in eq. (6). 

Fig. 6b Same as (a) at y/s = 500 GeV. 

Fig. 7 Parton polarizations Aq/q versus x taken at Q 2 ~ following 
eqs. (15, 16a, 16b). 

Fig. 8a The parity- violating helicity asymmetry Al versus y for W + and 

W~ production in pn collisions at y/s = 350 GeV. 
Fig. 8b Same as (a) at y/s = 500 GeV. 

Fig. 9a The ratio (see eq. (20)) versus y at = 350 GeV. Solid 
curve, 5q(x) =/ and Aq(x) =/ 0. Small dashed curve, 5q(x) ^ 
and Aq(x) = 0. 

Fig. 9b Same as (a) at ^ = 500 GeV. 

Fig. 10 The parity-violating helicity asymmetry A L versus y for Z pro- 
duction in pn collisions at y/s = 350 and 500 GeV. 

Fig. 11a The ratio af (see eq. (22)) versus y at y/s = 350 GeV. Solid 
curve, 5q(x) =/ and Aq(x) =/ 0. Small dashed curve, 5q(x) ^ 
and Aq(x) = 0. 

Fig. lib Same as (a) at y/s = 500 GeV. 
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Fig. 12a Parity-conserving double helicity asymmetry A LL versus y for W ± 

and Z production in pp collisions at y^s = 350 GeV. 
Fig. 12b Same as (a) at yfa = 500 GeV. 

Fig. 13 Double spin transverse asymmetry Att versus y for Z production 
in pp collisions at y/s = 350 and 500 GeV. 
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